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Abstract 

In the light of new experimental results on B — * Kir decays, we study the decay processes B — > Kit 
in the framework of both R-parity conserving (SUGRA) and R-parity violating supersymmetric 
models. We find that any possible deviations from the Standard Model indicated by the current 
data for the branching ratios and the direct CP asymmetries of B — ► Kir can be explained in both 
R-parity conserving SUGRA and R-parity violating SUSY models. However, there is a difference 
between the predictions of both models to the time-dependent CP asymmetry observable S^^o 
whose current experimental results include large uncertainties. We demonstrate that this difference 
can be useful for testing both models with more accurate data for Sk s7Z ° an d A^,p in the near future. 
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The quark level subprocesses for B — > Kit decays are b — > sqq (q = u, d) penguin processes 
which are potentially sensitive to any new physics effects beyond the Standard Model (SM). 
All the B — > Kit modes have already been observed in experiment and their CP-averaged 
branching ratios (BRs) have been measured within a few percent errors by the BaBar and 



Belle collaborations 
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The measurements of CP asymmetry observables 



r ere < 


3rr 


o 


rs 


so 


Q 


11 




12 


. ] 



decisive conclusions until recently |1|, |8j, |9j, 110|, llll LL2t LUl ■ But, the direct CP asymmetry in 
B° — > K ± tt t has been recently observed at the 5.7a level by BaBar and Belle 10, 11. 112I ] 
whose values are in good agreement with each other: the world average value is 

A^ P = -0.119 ±0.019 . (1) 

The direct CP asymmetry data for the other B — > Kir modes still involve large uncertainties: 
e.g., for B ± -> K^ modes, Agp = +0.04 ± 0.04. 

The recent experimental data for the CP-averaged BRs of B — > Ktx may indicate a 
possible deviation from the prediction of the SM: 

2B +0 „ B+- 



1.00 ±0.09 



R r , 



0.79 ±0.08 



(2) 



go+ ' n 2B 00 

where B 1 ^ denote the CP-averaged BRs of B — > K % ^ decays. It has been also claimed that 
w ittlin the SM, Rr Kfl „ flQ. The ahove expe^tal da t a s how the pattern H c > H. 
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11^ . which would indicate the enhancement of the electroweak (EW) penguin and/or 
the color-suppressed tree contributions 

On the other hand, in the conventional prediction of the SM, Aq° p is expected to be almost 
the same as Aqp'. in particular, they would have the same sign. However, the current data 
show that Aq P differs by 3.5a from A^p. This is a very interesting observation with the 
new measurements of A^p by BaBar and Belle, even though the measurements of Aq P still 
include sizable errors. This possible discrepancy from the SM prediction, together with the 
above one on R c and R n , has recently been called the U B — > Kit puzzle". One may need to 
explain on the theoretical basis how this feature can happen. 

In the light of those new data, including the direct CP asymmetry in B° — > K ± tt t : mauy 



works have been recently done to study the implications of the data 
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in id — > J\ ^7r^jjnany 



271 ] . However, most of those previous works have focused on finding 
out the B — > Kit puzzle itself and clarifying its implications through model-independent 
approaches, such as the topological quark diagram approach. 



factorization (QCDF) 



In this letter, we focus on how to resolve the B — > Ktt puzzle with well-motivated new 
physics models: in the framework of R-parity conserving and R-parity violating supersym- 
metry (SUSY). We calculate the BRs and the direct CP asymmetries for all the B — > Ktt 
modes in the SM and its SUSY extension with R-parity (SUGRA models) and without R- 
parity. Then, we present predictions of the different SUSY models to the mixing induced 
CP violating parameter Sk s tt° which has been observed with large errors through the time- 
dependent CP asymmetry measurement of B° — > K s tt° jil^J. In the recent work Q], it 
has been explicitly shown that the color-suppressed tree contribution is very sensitive to the 
observable Sk s -k q , while in contrast, the EW penguin contribution is not sensitive to S , ^ g7r o. 
As we shall see, the different SUSY models give different predictions to the time- dependent 
CP violating parameter Sir tt which can be tested by experiment. 

For calculation of the relevant hadronic matrix elements, we adopt the QCD improved 
28^ . This approach allows us to include the possible non-factorizable 
contributions, such as vertex corrections, penguin corrections, hard spectator scattering con- 
tributions, and weak annihilation contributions. The relevant end-point divergent integrals 
are parameterized as [3] 

*W. T s ( 1 + < 3 > 

where Xh and denote the hard spectator scattering contribution and the annihilation 

contribution, respectively. Here the phases <p HA are arbitrary, 0° < <p H A < 360°, p HA are 

free parameters to be of order one, typically p A < 2, and the scale A h = 0.5 GeV being the 

typical hadronic scale 2^ |. 

We first summarize the current status of the experimental results on B — > Ktt modes in 

Table HI which includes the BRs, the direct CP asymmetries (Acp), and the mixing- induced 

CP asymmetry (Sk^o). In order to exhibit the sign convention for CP asymmetries used 

in this work, let us specify the definition of CP asymmetries for B — > Ktt as follows. The 

direct CP asymmetry for B ± — > K ± tt° is defined as 

= B(B~ -> K~tt°) - B{B+ -> K+tt ) 
CP B(B~ -> K-7T ) + B(B+ -> K+7t°) ' 1 ' 

The definition of direct CP asymmetries for other B — > Kit modes becomes obvious. The 

time-dependent CP asymmetry for B° — > K s tt° is defined as 

= T(B°(t) -> K s *°) - T(B°(t) -> K s tt°) 
•* K s*°\ ' ~ T(B°(t) -> K s tt°) + T(B°(t) -> K s ir°) 
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TABLE I: Experimental data on the CP-averaged branching ratios {B in units of 10 -6 ), the direct 
CP asymmetries (Acp), and the mixing-induced CP asymmetry (S Ks7T o) for B — ► Kir modes. 
The SjCjrO is equal to sin(2^i]in the case that tree amplitudes are neglected for B° — > /T s 7r° 

,Q,y,y,Q,DJ3,Q- 



BR 


Average 


CP asymmetry 


Average 






24.1 ± 1.3 




-0.02 ±0.04 




- K±7r°) 


12.1 ±0.8 


Aq P 


±0.04 ± 0.04 




-> fl^Tr*) 


18.9 ±0.7 


A op 


-0.115 ±0.018 


B(B° 


-> ^ vr ) 


11.5 ± 1.0 




±0.001 ±0.155 
±0.34 ±0.29 



= S K n o sin(Am d t) - C K n o cos(Am d t) 



(5) 



where T denotes the relevant decay rate and Amj is the mass difference between the two 
B° mass eigenstates. The S Ks7r o and C Ks7r o are CP violating parameters. In the case 
that the tree contributions are neglected for B° — > K S 7T°, the mixing-induced CP violating 
parameter Sk s k° is equal to sin(20i) [fa (= j3) is the angle of the unitarity triangle]. Note 
that the measured value of w o (Table |l} is different from the well-established value of 
sin(20i) = 0.725 ± 0.037 measured through B -> J/ipK^ 1]. It may indicate that the EW 
penguin and the color-suppressed tree effects play an important role jlfl ]. 

In the following two sections, we will discuss possible resolutions of the B — > Kit puzzles 
in the context of SUSY models. 



[1] R-parity violating SUSY case 

In the R-parity violating (RPV) minimal supersymmetric standard model, we will assume 
only /'—type couplings to be present The R-parity violating interaction introduces new 
operators. The relevant new operators are 

Ceff = (Ual^Up) (SffypRba) - W^Wi) ^L(R)d p ) (d^R(L)b a ) 

(6) 

(LlMm) (WpR(L)b a ) , (7) 
where L(R) = (1 ±7 5 )/2, a and (3 are the color indices, and mj denotes the sfermion mass. 
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Note that the operators having the following chirality structure (paj^Lqp) {fp^^Rba) do not 
exist in the SM effective Hamiltonian. 

The RPV SUSY part of the decay amplitudes of B — > Ktt modes are given by j^] 



A RPV (B° - 
A RPV (B- 



> K-TT°) 



-if K F B ^(0) (m 2 B - ml) uf 12 R K c A + A^{K~^ 
iUF B ^ K (0) (ml -ml 



(8) 



Uu2 ^/2 (~ TK7tRkCa + a ') 



^U2 ^121) ^^Rtt c a (df 2 i ^112) / 



.4 



RPV 



B - -> K\-) = zf K F B ^(0) (m\ 



m„ 



V2 



(diu ~ ^121) a ' - (d-ui - ^112) RkCa 



(9) 
(10) 



A RP V( £0 _ R o n0) = iUF *-+K(p) ( m l _ 



K 



X 



V2 
(d R -d L 



V2" 



-r Kn R K CA + a) 



where the annihilation contributions are given by 



7T , 



= -iIbUk [(d&j, - d L ni ) b' 4 + (df 21 - 4 12 ) &' 3 ] (12) 
-ifsUfK uf 12 b' 3 . (13) 



The u?u„ and cL»lf are defined as u^i 



l jkn 



jkn 



jkn 



V J 1 * 

E3 l ijn\k-A JR 
i=l 8^1 ' "ifc" 



E3 ^iZk^inj 
i=l 8m? 



We refer to Refs. [2j| for the relevant notations. Here fi and F B ~* 1 denote decay constants 
and form factors, respectively. The parameters a', Ri, are defined as 



a' = 


ca 


[1 Cf(Xs v' p 


ca 


Rk 




2m 2 K 




m h (p)(m q (p) + m s (/i)) 


R n 




2ml 

TV 




m b (n)(m u (fi) +m d (/i)) 



s w 

n PiPx ' 



u (d) for K~ (K )) 



(14) 
(15) 
(16) 

UF B ^(0)(ml-m%) ' (17) 
where iV c (= 3) is the number of colors and Cp = (N% — 1)/(2N C ). Vp 2 and H' PiP2 come 
from the vertex corrections and the hard spectator scattering contributions, respectively. 
For their explicit expressions, we refer to j3l|. P\ is the final state meson absorbing the 



f K F B ^(0)(m%-mi: 
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FIG. 1: R n versus R c (left one) and A~^, p versus A^° p (right one) in the R-parity violating SUSY 
model. 

light spectator quark from B meson and P 2 is the other final state meson emitted without 
absorbing the spectator quark. The parameters b\ are defined as 

&3 = % c c c 4 , K = % c C c A{ , (18) 

where 

2 \ 



A\ = ira s 



18 (x A -A + ^-\+2r 2 x X 2 A 



A{ = \2ita s r x (2X\-X A ) , (19) 

with r x « R v . X A is the divergent integral as defined in Eq. (jSJ). c A c are the RGE improved 
QCD enhanced factors at the scale fi = m^. 

From Eqs. © — (fTT| . we note that the R-parity violating couplings df- k and d\ mn always 
appear as the combinations (df 12 — d\ 2X J and (df 21 — df 12 ). Thus, in this analysis, we ac- 
tually use three different combinations of R-parity violating couplings: wf 12 , (^m — ^m) 
and (c?f 21 — rff 12 )- Since each combination can be expressed as a complex number, we have 
six independent real parameters arising from the new physics effects and we have 9 results 
to explain. The contributions of the new terms to the amplitudes are mostly different for 
different decay modes. 

By varying the above parameters, we try to fit all the current data simultaneously as 
shown in Table |TJ In Fig. [H we show R n versus R c (left figure) and A^p versus Aq° p (right 




R c -R n A cp [KV] 

FIG. 2: Sk n o versus (R c — R n ) (left one) and Sj^ n o versus A^,p (right one) in the R-parity 
violating SUSY model. 

figure). Here the same parameter sets are used to fit both the BRs and the direct CP 
asymmetries. We see that the values of R n , R c , A^p, and A^° p are consistent with the 
current data at la level. In fact, it turns out that all the current data for the BRs and the 
direct CP asymmetries, including A^p and A^p, can be explained at la level in the R-parity 
violating SUSY model. In other words, the possible discrepancy between the SM predictions 
and the current data for the BRs and the direct CP asymmetries can be explained by the 
new physics contributions which, in particular, come from the new operators having the new 
chirality structure as mentioned below Eq. Q. 

Using the same values of the parameters used in Fig. ^ we predict the mixing induced 
CP violating observable Sk s tt°- In Fig. El our result is presented as Sk s tt° versus (R c — R n ) 
(left figure) and Sk s -k° versus A^p (right figure). We see that our prediction is in good 
agreement with the current data at la level. We shall see in next section that in R-parity 
conserving SUSY case, it is very difficult to explain the small value of the current data for 
S Ks7r o together with the other data, especially R c and R n . 

In Table m we show the representative values of our prediction to the BRs, the direct CP 
asymmetries and the mixing induced CP asymmetry in the R-parity violating SUSY model. 
We consider two cases: (i) p HA = and (ii) p A = 0.3, p H = 4> HA = 0. The corresponding 
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TABLE II: Predictions of the R-parity violating SUSY model for two cases: (i) p H A = and (ii) 
p A = 0.3, p H = (pH,A = 0. The case (ii) are shown in the bracket. (B in units of 10~ 6 ) 



BR 


Prediction 


CP asymmetry 


Prediction 


B(B ± -> J RT°7r ± ) 
B(B ± i^vr ) 
S(S° -► K^iF) 

B(B° -» ^°7T ) 


23.6 [24.8] 
13.3 [13.0] 
19.0 [18.9] 
11.9 [12.6] 


j\ CP 
A+0 

aOO 

J\ CP 

Sk s tt° 


-0.010 [-0.007] 
+0.026 [+0.018] 
-0.134 [-0.115] 
-0.142 [-0.141] 
+0.51 [+0.55] 



values of the couplings are (in 10 8 ) 

Ki 2 -4 J1 |~3.2 (3.1), |df 21 -df 12 |~ 0.87 (0.60), |< 12 | ~ 2.1 (2.2) (20) 

The values in the parenthesis are for the p A = 0.3 case. The constraints on the RPV 
couplings need to be checked. However, apart from wf 12 , the rest of the couplings appears in 
the amplitude as combinations (e.g., df 12 — d^ 2 i) of 3 or 4 different RPV couplings A- jfc so that 
they easily satisfy the constraints. uf 12 involves A- 12 A^ 3 . In our example above (for p H A =0 
case), A^ ~ 8 x 10 -2 was used. It is also important to note that uf 12 involves mf which we 
assume to be ~ 200 GeV. The experimental bound on X' 3lk is given by A' 3U , < 1.2 x 10 _1 for 
1 TeV of squark mass by using the ratio of BRs of K + — ► 7i + uu and K + —> 7r°z/e + decay j^]. 
However, the bound on A' determined from the experimental value of the BR of K — > -nvv 
decay depends on the squark mass and in GUT models, it is quite natural to expect a large 
hierarchy (~ 5) between the squark and the slepton masses. 
[2] R-parity conserving SUSY case 

In this case the SUSY contributions appear in loop. The one loop SUSY contributions 
are available in the literature, e.g., Refs. [yj]. In our calculation, we do not use the 
mass insertion approximation, but rather do a complete calculation. The SUGRA model 
starts at the GUT scale. We assume the breakdown of the universality to accommodate the 
B — > 71K data. While we satisfy this data, we also have to be careful to also satisfy other 
data, e.g., b — > 57, AMk, ABj, €k, etc. 

We use the following boundary conditions at the GUT scale: 

( m % LL ,u RR , DRR) = ml (6» + ^ ll , Urr , Drr) ) ; 4^ = A (Yfa + A^ d) ) . (21) 
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TABLE III: Predictions of the R-parity conserving SUSY model. The SUSY parameters are men- 
tioned in the text. (B in units of 10 -6 ) 



DiX 


Prediction 


CP asymmetry 


Prediction 


B(B ± -» K ^) 
B(B ± -» i^vr ) 
B(B° -> #±7^) 
£(5° £T 7t°) 


23 
10.3 
19.1 
11.3 


40+ 

"^CP 
•^CP 

A 00 

Asp 


-0.030 
-0.0073 
-0.105 
-0.08 
+0.73 



The SUSY parameters can have phases at the GUT scale: m» = |mi/ 2 |e 1 ^ (z = 1,2,3) 
(the gaugino masses for the U(l), SU(2) and SU(3) groups), A = \Ao\e ta A and \i = 
|/z|e ?6,M . However, we can set one of the gaugino phases to zero and we choose 82 = 0. The 
electric dipole moments (EDMs) of the electron and neutron can now allow the existence 



of large phases in the theory 



35 



36, 371. In es but 

n 

calculate the EDMs to make sure that current bounds (\d e \ < 1.2 x 10 27 ecm [38] and 
\d n \ < 6.3 x 10~ 26 ecm Q) are satisfied. 

We evaluate the squark masses and mixings at the weak scale by using the above boundary 
conditions at the GUT scale. The RGE evolution mixes the non-universality of type LR (A 
terms) via drriQ 2 LLRR /dt oc A^ u ^A u ^) terms and creates new LL and RR contributions at 
the weak scale. We then evaluate the Wilson coefficients from all these new contributions. 
We have both chargino and gluino contributions arising due to the LL, LR, RL, RR up type 
and down type squark mixing. These contributions affect the following Wilson coffecients 
C3 — C10, C7 7 and C$ g . The chargino contributions affect mostly the electroweak penguins 
(C7 and Cg) and the dipole penguins, while the gluino penguin has a large contribution to the 
dipole terms due to the presence of an enhancement factor rrig / rrib (the gluino contribution 
also affects the QCD penguins, but the effect is small). We include all contributions in our 
calculation. The SUSY contributions also bring new operator contributions over the SM by 
having a chirality exchange in the SM operators. 

The electroweak penguin contribution is required to solve the B — > nK puzzle for the 
BRs and can solve the CP asymmetries Q]. If we do not consider the BRs, then the direct 
CP asymmetries of the B — > ttK modes can be solved by the dipole penguin contributions 
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-0.1 

a cp [kV] 
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FIG. 3: S K ^o versus (R c — R n ) (left one) and SV^o versus A^, p (right one) in the SUGRA model. 

only. The dipole penguin contributions can not be arbitrarily large, since it is also present 
in the b — > 37. In order to obtain a fit, we find that A^ d are necessary. The nonzero values 
of these parameters generate the dipole penguin and the (Z-mediated) electroweak penguin 
diagrams. In Table ITTT1 we show an example of a fit. From the fit one finds the prediction for 
S Ka ^o to be large. The SUSY parameters used for this fit are: m 1( / 2 = 450 GeV, A = —800 



GeV, m Q = 300 GeV, A 2 . 3 = 0.2 e" - 3 *, AAf = 0.55 e a8 \ AAf = 0.05 e~ L5 \ tan/5 = 40, 



H > 0. Since the SUSY parameters have phases, the EDMs of the electron and the neutron 
need to be checked, and we do indeed satisfy the experimental bounds for these EDMs. For 
this example, we find \d e \ = 2.23 x 10~ 29 e cm and \d n \ = 8.2 x 10 -27 e cm. The QCD 
parameters for this fit are: p A = 2 and (ft A = 2.77. In this fit we have used nonzero Aq ll , 
but it is possible to obtain fits without Aq ll . We can obtain fits for other tan (3 values as 
well. 

In Fig. 3, we show Sk s k° as a function of (R c — R n ) and Sk^q as a function of A^p. In 
order to generate these figures, we have varied mi/2, mo, tan/5 and A's. We see from the 
figures that the lowest value of Sk s7 t° is about 0.69 and the maximum direct CP asymmetry 
Agp predicted by the SUGRA model is about —0.107. If we compare Figure 3 with Figure 
2, we find that the prediction for Sx g7r o in the R-parity conserving SUSY model is much 
higher than in the R-parity violating SUSY model and therefore the future data on SV s7r o 
will be crucial. The future data (with reduced error) of A^p is also crucial to distinguish two 
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scenarios since the maximum direct CP asymmetry A^p predicted by the SUGRA model 
is about —0.107, whereas the asymmetry can be larger negative in the R-parity violating 
model. 

In conclusion, we have explained the recent experimental results on the BRs and CP 
asymmetries of different B — > ttK modes in R-parity violating and R-parity conserving 
SUSY models. We have found that the R-parity conserving SUSY model tends to generate 
large Sk s tt° when we use all the constraints on the BRs and CP asymmetries, and the lowest 
value of Sk s tt° is about 0.69. However, lower values of Sk s -k° can be accommodated in the 
R-parity violating SUSY model. We also find that the future data of A^p is important to 
distinguish the two models. 
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